PAIRING MECHANISM IN THE TWO DIMENSIONAL HUBBARD MODEL 
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Here we discuss Quantum Monte Carlo results for the magnetic susceptibility, single— particle spectral weight 
and the irreducible particle-particle interaction vertex of the two-dimensional Hubbard model. In the doped 
system, as the temperature is lowered below J = 4i 2 /U , short— range antiferromagnetic correlations develop. 
These lead to a narrow low-energy quasiparticle band with a large Fermi surface and a particle-particle 
vertex which increases at large momentum transfer, which favor d I 2_„2-pairing. 



g : A ™iety of experiments^ the high T c cuprates can be interpreted in terms of a gap with 



dominant d x 2_ y 2 symmetry |]J]. However, the implications of this with respect to the nature of 
the underlying pairing mechanism remains an open question. Here we review some results which 
have been obtained for the two-dimensional Hubbard model, which has been found to exhibit 
d a .2_j ( 2-like pairing fluctuations. Our aim is to examine in this particular case the structure of 
the quasiparticle spectrum and the interaction in order to gain insight into the mechanism which 
leads to d x 2_ y 2 pairing correlations in this model. 

The two-dimensional Hubbard Hamiltonian provides an approximate model of a CuC>2 layer: 



2 ■ H = -t ]T {c\ s c j8 + cJ.Ci.) +UJ2 nw. (1) 

^ ! «>,» i 

Here c is creates an electron of spin s on site i, and n, s = c is Ci S is the occupation number for 
spin s on site i. The one-electron transfer between near-neighbor sites is t, and U is an onsitc 
Coulomb energy. The bare energy scale is set by the bandwidth 8t and the effective Coulomb 
interaction U, which are both of order electron volts. Near half-filling, electrons on neighboring 
sites tend to align antifcrromagnetically so as to lower their energy by the exchange interaction 
J = 4t 2 /U. This interaction is of order a tenth of an electron volt and, as we will see, sets the 
, energy, or temperature scale, below which antiferromagnetic (AF) correlations, the low-energy 

structure in the single-particle spectral weight, and the pairing interaction develop. 

While Monte Carlo S and Lanczos J^] calculations for a 4 x 4 lattice find that two holes 
added to the half-filled Hubbard ground state form a d x 2_ y 2 bound state, and density matrix 
renormalization group calculations Q find that d x 2_ y 2-\ike pairs are formed on two-leg Hubbard 
ladders, it is not known what happens for the two-dimensional Hubbard model. It is possible 
that on an energy scale of order J/10, a d x 2_ y 2 superconducting state forms. However, this may 
well require modifications of the model, such as an additional near-neighbor A JSi ■ Sj term or 
possibly a next-near-neighbor hopping t'. Nevertheless, it is known that as the temperature 
is reduced below J, d x 2_ y 2 pairing correlations develop in the doped two-dimensional Hubbard 
model, and here we will examine why this happens. 

At half-filling, (n^ +71^) = 1, the 2D Hubbard model develops long-range antiferromagnetic 
order as the temperature goes to zero. In the doped case, strong short-range AF correlations 
develop as the temperature decreases below J. This is clearly seen in the temperature dependence 
of the wave vector dependent magnetic susceptibility 

X(<0 = f dr (mr + ,(rK+(0)>e-^. (2) 
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Figure 1: (a) Magnetic susceptibility x( ( l) versus q along the (1,1) direction for various temperatures, 
(b) Temperature dependence of the AF correlation length £af in units of the lattice spacing a. These 
results are for an 8 x 8 lattice with U/t = 4 and a filling (n) = 0.875. 



Here mf = ct-.Cj| and Jn7^_^(r) = e HT m^ +p e~ HT , where m~ +e is the hermitian conjugate of m^ +p . 
Monte Carlo results for xil) versus q along the (1, 1) axis for an 8 x 8 lattice with U/t = 4 and 
a filling (n) = 0.875 are shown in Fig. [l](a). As the temperature decreases below J = 4t 2 /U, 
significant short-range dynamic AF correlations evolve. Fig. |l](b) shows the AF correlation 
length £af versus T. Here, is denned as the half-width at half-maximum of x(q). 
As these AF correlations develop, the single-particle spectral weight || 

A(p, u>) = Im G(p, iuj n — > lo + iS) (3) 

and the density of states N(ui) — ^ p A(p, u>) also change. Figure ||(a) shows N(u>) for U/t = 8 
and (n) = 0.875. As the temperature is lowered, a peak appears on the upper edge of the lower 
Hubbard band. 

This peak arises from a narrow quasiparticle band shown in the single-particle spectral weight 
A(p,oj) of Fig. H and plotted as the solid curve in Fig. ||(b). As the momentum p goes towards 
the r point (0,0), we believe that the quasiparticle peak is obscured by the lower Hubbard 
band because of the resolution of the maximum entropy technique which we have used. Indeed, 
at the r point a separate quasiparticle peak is found from Lanczos exact diagonalization on a 
4x4 lattice Q. Note that this implies a large hole-like Fermi surface. As clearly evident in 
the spectral weight shown in Figs. ^|(a) and (b), the quasiparticle dispersion is anomalously flat 
near the (it, 0) corner. As discussed by various authors, this reflects the influence of the AF 
correlations on the quasiparticle excitation energy. It is clear that the peak structure in N(to) 
also arises from the short-range AF correlations and is a many-body effect rather than simply a 
non-interacting band Van Hove singularity. 

Monte Carlo calculations have also been used to determine the singlet irreducible particle- 
particle vertex rig(p', —p',p, —p) in the zero center-of-mass momentum and energy channel which 
gives the effective pairing interaction. Here p = (p, iu) n ). In Fig. 4(a), Tig(q — p — p') is plotted 
for q along the (1, 1) direction and iuj n — iuj n r — iirT, corresponding to zero Matsubara energy 
transfer. Comparing Figs. 1(a) and 4(a), one clearly sees that the structure of the interaction 
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Figure 2: (a) Evolution of the single-particle density of states with temperature for U/t = 8 and 
(n) = 0.875. (b) Dispersion of the quasiparticle peak in the spectral weight versus p. The solid points 
mark the low-energy peaks of A(p,io) shown in Fig. 3, and the solid curve represents an estimate of the 
quasiparticle dispersion using these data and Lanczos results for p near (0,0). The broad darkened areas 
represent the incoherent spectral weight in the upper and lower Hubbard bands. The horizontal dashed 
line denotes the chemical potential fx. 
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Figure 4: (a) Singlet irreducible particle-particle vertex for zero energy transfer Fis (q, — 0) versus q 
along the (1,1) direction. As the temperature decreases below J = 4t 2 /U, the strength of the interaction 
is enhanced at large momentum transfer. Note the similarity to x(q) m Fig- l( a )- 3 ) Temperature 
dependence of the d x 2_ y 2 eigenvalue of the Bethe-Salpeter equation. These results are for U/t — 4 and 
(n) = 0.875. 



and x(q) are similar, both reflecting the development of the AF correlations as T is reduced 
below J. 

Given the Monte Carlo results for the irreducible particle-particle vertex Tjsip', —p',p, —p) in 
the zero energy and center-of-mass momentum channel, and the single-particle Green's function 
G(p,iui n ), the Bethe-Salpeter equation for the particle-particle channel is 

\Mp) = -^^r IS (p,-p,p',-p')|G(p')l 2 ^(pO- (4) 

Here, the sum on p' is over both p' and uj n >. In the parameter regime that the Monte Carlo 
simulations are carried out, the leading singlet eigenvalue is in the d x 2_ y 2 channel. Fig. 4(b) 
shows the temperature dependence of the d x 2_ y 2 eigenvalue. 

The development of both the low-energy quasiparticle dispersion and the peak in the singlet 
particle-particle vertex at large momentum transfers arises from the growth of short-range AF 
correlations as T decreases below J. As is known [pj, for the large Fermi surface associated 
with the observed quasiparticle dispersion, a particle-particle vertex which increases at large 
momentum transfer favors d x 2_ y 2 pairing. Note that the tendency for d x 2_ y 2 pairing does not 
require a particularly sharp or narrow peak in ris(q) for q = (jr, 7r), but rather simply sufficient 
weight at large momentum transfers. Thus it is the strong short-range AF correlations which 
lead to the formation of d x 2_ y 2 pairing correlations in the Hubbard model. 
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